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We propose an anharmonic oscillator driven by two periodic forces of different frequencies as a 
new time-dependent model for investigating quantum dissipative chaos. Our analysis is done in the 
frame of statistical ensemble of quantum trajectories in quantum state diffusion approach. Quantum 
dynamical manifestation of chaotic behavior, including the emergence of chaos, properties of strange 
attractors, and quantum entanglement are studied by numerical simulation of ensemble averaged 
Wigner function and von Neumann entropy. 
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Ph ! I. INTRODUCTION 

_^ Quantum nonlinear systems with chaotic classical counterparts have received much attention in the last two decades. 

I This field of investigation is sometimes called "Quantum Chaos" 0. The usual procedure of studying a quantum 
chaos is to take a system which exhibits chaotic motion, when treated classically, and see what effects occur in a 
I ' quantum treatment. All real quantum systems are open and their classical limit is related to the loss of coherence 
^ ' produced by interaction with the environment. Thus, investigations of quantum chaotic system are connected with 
\l ■ the correspondence problem in general, and with decoherence and dissipation in particular. 

] ] It is now well established that the quantum dynamics of classically chaotic systems will show major departures from 
. the classical motion on a suitable time scale. Among these phenomena we note the dynamical localization of classical 
■ diffusive excitation, due to quantum mechanical interference, which is in close analogy with the Anderson localization 
in a random potential. Dynamical localization has been well studied theoretically |^ and verified in experiments |^ 
with laser-cooled atom moving in a standing wave with periodically modulated nodal position. Much research on the 
subject of classical and quantum chaos is devoted to the kicked rotator, which exhibits regions of regular and chaotic 
motion in the Poincare section (see for example [Q). This model is very popular in investigations of transition to 
quantum chaos. Its experimental realization, and observation of the model's dissipation and decoherence effects, is 
^ . carried out on gas of ultracold atoms in a magneto-optical trap subjected to a pulsed standing wave In Ref. 

d ' it was proposed to realize the parametrically kicked nonlinear oscillator model in a cavity involving Kerr nonlinearity. 
^ , It was also shown that more promising realization of this system also including quantum regime, is achieved in the 
dynamics of cooled and trapped ion, interacting with periodic sequence of both standing wave pulses and Gaussian 
^ ^ laser pulses Another suggestion to investigate quantum chaos in a single trapped ion was recently provided in 
1^. In fact, while numerous theoretical works on the subject of quantum chaos were carried out, their experimental 
realization remains somewhat scarce. Moreover, there are still many questions to be answered, yet, and there is 
. obvious need for new systems showing chaos, and new experiments. 
" " ' In this paper we propose a new kind of physical systems showing dissipative chaotic dynamics. These systems are 
modeled by a dissipative nonlinear oscillator driven by two periodic forces of different frequencies. This model was 
proposed to study quantum stochastic resonance in the authors' previous paper |l0| , where it was shown, in particular, 
that the model is available for experiments. It can be implemented at least for dynamics of strongly interacting photons 
in optical cavity with x*-^^ nonlinear medium ||Tl[| , and for cyclotron oscillations of a single electron in a Penning trap 
1^,^. We would especially prefer to study the transition to quantum chaos and its control, the role of dissipation 
and quantum entanglement in chaotic dynamics, and adjacent questions of characterizing quantum chaos. These 
investigations are complemented by consideration of the information aspect of chaotic dynamics through the von 
Neumann entropy. 

In classical mechanics standard characterization of chaos might be given in terms of the unpredictability of phase- 
space trajectories. However, the most important characteristics of classical chaotic systems - exponential divergence 
of trajectories, starting at arbitrarily close initial points in phase-space - does not have quantum counterpart. The 
question of what the quantum mechanical equivalent of chaos constitutes, has been posed. Many criteria have been 
suggested to define chaos in quantum systems varying in their emphasis and domain of application As yet, 

there is no universally accepted definition of quantum chaos. Our analysis of quantum chaos is based mainly on the 
time-evolution of von Neumann entropy and Wigner function. 
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The system under research is dissipative and therefore it has a mesoscopic nature. Quite generally, chaos in classical 
conservative and dissipative systems has completely different properties, e.g. strange attractors can appear only in 
dissipative systems. Therefore, the system of our interest might allow us to examine challenging problems of quantum 
dissipative chaos, including the problem of quantum counterpart of a strange attractor. We note, while the quantum 
dynamics of isolated or so called Hamiltonian systems with chaotic classical counterparts has studied well, a very 
liitle work has been done in looking at quantum chaos of dissipative nonlineat systems. Among earlier studies of 
open quantum chaotic systems it should be noted the papers [l5| ]. A new impetus to increasing the dissipative 
quantum chaos area has recently given by study of decoherence and quantum-classical correspondence problem of 
chaotic systems p6| , p7[ . 

Below we use the traditional ensemble description of Markovian open systems, based on the master equation. Then, 
this equation is presented in quantum trajectories in the frames of quantum state diffusion approach (QSD) [ISj. 
Recently, it was shown how quantum state diffusion can be used to model dissipative chaotic systems on individual 
quantum trajectories p9| , p0[ |. In contrast with these papers, here we show how it is possible to describe the quantum 
chaos using a statistical ensemble of trajectories, which is actually realized in nature. Our results indicate, that 
though properties of quantum chaotic dynamics do not appear on ensemble averaged oscillatory excitation numbers, 
these properties are clearly obvious on the entropy and Wigner function. 

The outline of this paper is as follows: In the next section we describe the system proposed, and give the analysis 
of its classical motion on the Poincare section in phase-space. In Sec. Ill we develop the quantum description of the 
problem by numerically solving the master equation through the QSD method. We present the averaged over quantum 
trajectories results for the mean excitation number of nonlinear oscillator, and for both the Wigner function and the 
von Neumann entropy. We summarize our results in Section IV. 



II. MODEL OF DOUBLE DRIVEN OSCILLATOR: CLASSICAL PHASE MAP 



In this section we give the theoretical description of the system. The nonlinear oscillator driven by two periodic 
forces at frequencies ui and 102 and interacting with a reservoir is described by the following Hamiltonian 

H — hujQa'^a + ftx(a^a)^ -|- h [{fli exp(— iwit) -|- exp{—iu}2t)) + /i.e.] -|- Hioss, (1) 

where a, a"*" are boson annihilation and creation operators, is an oscillatory frequency and x is the strength of the 
anharmonicity. The couplings with two driving forces are given by Rabi frequencies Oi and O2. -ffjoss = oT^ + a~''r 
responsible for the linear losses of oscillatory states, due to couplings with heat reservoir operators giving rise to the 
damping rate 7. The reduced density operator p within the framework of the rotating- wave approximation, in the 
interaction picture corresponding to the transformation p — » e~*"^° at p^iuiia at governed by the master equation 

% = -\ + H^ntA + {^^P^t - \LtL,P - \pLtL^ , (2) 

where 

Ho = hAa+a, (3) 
Hint = ft [{^1 + ^2 exp {-iSt)) a+ + (fij 4- exp (iSt)) a] + hx{a^af. 

Here A. — loq— uji is the detuning, and 5 = u}2— i^i is the difference between driving frequencies, which works 
as modulation frequency. Li are the Lindblad operators: 

Li = ^{N + l)-fa, L2 = ^/N^a+, (4) 

where 7 is the spontaneous decay rate of the dissipation process, and N denotes the mean number of quanta of a heat 
bath. Note, only the case of a vacuum reservoir, iV = 0, considered below. 

For ri2 — Q this equation describes the single driven, dissipative anharmonic oscillator, which is a well-known and 
archetypal model in nonlinear physics |^,^. In the semiclassical approach and in steady-state this system exhibits 
bistability, which appears as the hysteresis behavior of the mean oscillatory number, versus either the detuning 
A or the strength of driving Q,i However, the hysteresis in the quantum-mechanical treatment disappears on 

ensemble-averaged mean oscillatory number n{t) = {a^ [t)a{t)) , and bistability, manifests itself on individual quantum 
trajectories, as noise-induced transitions between two possible metastable states p3| , as well as on the statistics of 
oscillatory numbers p^ . 
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In the case of double driven oscillator, when two external forces are both present, the corresponding Hamiltonian 
includes explicit time-dependence, even in the rotating-wave approximation. It may therefore, be expected that 
the system presented above, exhibits regions of regular and chaotic motion depending from the parameters: Xi A, fii, 
il2 and 7. To illustrate the operational regimes of the oscillator, first we pay attention to the classical description. 
After making the usual approximations, the classical limit of Eq.(^ becomes 

= -^7^ - i (a + x(l + 2 |a|')) a - 2 (fii + cxp {-iSt)) , (5) 

where a is the dimensionless complex amplitude, corresponding to the operator a. We analyze the time-dependent 
solution of this equation in the phase-space of dimensionless position and momentum X = Re a, Y = Ima. There 
are a number of ways to make such analysis. We adopt a discrete surface or Poincare section of this system. Let Xq, 
Yq be an arbitrary initial phase-space point of the system at the time to- Then we define a constant phase map in 
the {X, Y) plane by the sequence of points {Xn, Yn) = {X{tn), Y{tn)) , changing the time intervals by tn — to + 
n = 0, 1, 2, ... . This means, that for any t = tn the system is in one of the points of Poincare section. Our analysis 
shows that for extended time scales exceeding the damping rate, the asymptotic dynamics of the system is regular 
in the limits of small and large values of modulation frequency, i.e. ^ <C 7 and 5^7, and also when one of the 
perturbation forces is much greater than the other: fli ^ or ^2 ^i- The dynamics of the system is chaotic in 
the range of parameters S > j and fii cf.il.2- Fig.l shows the results of numerical calculations of the classical maps, 
for the parameters chosen in the range of chaos. These figures clearly indicate the classical strange attractors with 
fractal structure, which are typical for a chaotic Poincare section (we choose to = for all cases). It is expected from 
these results, that the domain of phase-space which includes an attractor, strongly depends on the parameters fii, 
and 7. In particular, the domain increases with the increase of fli or 1^2, and with the decrease of the decay rate 
7. Below we give the general consideration of this effect considering equation (pf) in the following integral form: 



a{t) = 



(6) 



-i j {ni + n2e~'^'')e'f^''^e^*'/^dti+a{to)e^*"/' 

to 

Here ao — a(to) is an initial value at t — to, and f{t) is introduced as follows : 

t 

f{t)^ j (A + x(l + 2|a(ti)|^))dii. (7) 

to 

It is easy to estimate the module of the complex amplitude \a{t) \ , taking into account that f{t) is a real function. 
For its maximum on time interval value amax = maxf \a{t)\ we obtain 

^ (l^l| + l^2|) 

7 

This formula determines the border of a classical map in dependence from Rabi frequencies and the decay rate, and 
explains the above conclusion about the size of the strange attractor in phase-space. It is interesting, that this border 
size for the anharmonic oscillator is independent from the parameter of nonlinearity x and the detunings A and 5. 
It should be noted a peculiarity of the system proposed regarding the strange attractors. As can be seen from the 
numerical results, the attractors depicted on Fig. 1 (a) and Fig. 1 (c) for different parameters have the same form 
in phase-space and differ from each other only on the scales. It is obvious that this property is the consequence of a 
scaling symmetry of the classical equation ^ . Indeed, it easy to verify that Eq. (||) remains invariable for the following 
scaling transformation of the complex amplitude a — > a' = Aa, where A is a real positive dimensionless coefficient, 
if the parameters x, A, ^2 correspondingly transformed as: x ^ x' = X / ^ ^ ^' = ^ + X (l ^ 1 / '^^) : 
fii 2 — > 2 = ^^1,2- We have illustrated such symmetry on Figs. 1 (a),(c) where the attractors are presented for two 
sets of parameters coupled by the scaling transformations. We note, that in the quantum treatment the diffusion term 
in the master equation affects the scaling symmetry roughly. In Sec. Ill, we will discuss this and other properties 
of strange attractors common with the Wigner function, in more detail. 



III. QUANTUM SIGNATURES OF CHAOS: ENTROPY AND WIGNER FUNCTION 

In this section we examine the problem of quantum dissipative chaos on the basis of the von Neumann entropy and 
the Wigner function. The peculiarities and advantages of such a consideration for interpreting quantum chaos in an 
ensemble theory, are as follows: 
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The von Neumann entropy, which is defined through the reduced density operator as 



S=-Tr[p\np), (9) 

is a measure of dissipation and decoherence. The entropy for an isolated quantum system, does not change under 
the time evolution. If the time evolution of the system is perturbed through interaction with an environment, the 
averaging over the perturbation typically leads to an entropy increase. The von Neumann entropy also is a sensitive 
operational measure of an entanglement, as well as the measure of the purity of quantum states If the system 
is in pure state, the entropy is precisely zero. Thus, it is expected that the study of entropy for doubly driven 
oscillator will allow to examine chaos in terms of quantum entanglement. 

For an extended time-evolution period the classically chaotic system at moments <„ = to + "^n, rt = 0, 1, 2, fills 
the Poincare section in phase-space. There are a number of methods to treat Poincare section quantum mechanically. 
We note the Ref. pffl, where a method of quantization of classical dissipative maps was proposed. The Poincare 



section of one quantum trajectory was considered in [|19 20|. For our dynamical model we adopt the method of the 



Wigner function and study the features of correspondence between the Wigner function, which is taken at one of the 
moments t — tn — tt^ + and Poincare section. It was obtained through a computer simulation of the Wigner 
function. 

We analyze the problem of dissipation on the basis of the QSD method, which operates with stochastic states 
|\E'j(t)) , describing the evolution along a quantum trajectory. The equation of motion is: 

= {Ho + H,^t) I*?) dt - (10) 
\ J2 {LtU - 2 (£+) + (L,) 1*^) dt + Y. {U - {U)) 

i i 

Here ^ indicates the dependence on the stochastic process, the complex Wiener variables dS^i satisfy the fundamental 
correlation properties: 

M (de.) = 0, M {didij) = 0, M {didCj) - 5,jdt, (11) 

and the expectation value equals {Li) = {'i> ^ \Li\ ^) . According to this method the reduced density operator is 
calculated as the ensemble mean 

-J m 

p{t)=M{\^.){^^\)= hm -Y^\^^{t)){^^{t)\ (12) 

over the stochastic pure states |5'^(t)) , describing the evolution along a quantum trajectory. 

Let us first qualitatively describe the most important physical processes, which are responsible for the origin of 
chaotic dynamics in our QSD numerical study. There are two ways to realize the controlling transition from the 
regular to chaotic dynamics, by changing one of the parameters of the system. One of these ways is to vary the 
strength of the second force in the range from ^ to ^ ^i. In the limit ^ ^^i the system is reduced 
effectively to the model of single driven anharmonic oscillator, which exhibits bistability for the definite range of 
parameters x. A, f2i, and 7. In this limit our analysis of the time-dependent stochastic trajectories for expectation 
numbers n^{t) = (^'^ |a+a| ^'j) shows, that the system in the bistability range spends most of the time close to one 
of the semiclassical solutions of Eq.(|^), with quantum interstate transitions, occurring at random intervals. With 
increasing the amplitude Q2 for the other parameters, chosen to lead bistability in terms of semiclassical solution, 
the stimulated processes, i.e. dynamical interstate transitions, become sufficient. Their contribution at Q2 — ^1 
leads to the emergence of a chaotic regime on quantum trajectories. In the limit fl2 3> r2i the regular dynamics is 
restored. This limit is equivalent to the case of Qi ^ Q2, because we can choose such an interaction picture, that 
time dependent exponent in (^) will appear near Qi. In fact, using the transformation p g-«i^2a+at^gicj2a at Qf ^j^g 
reduced density operator, we arrive to the master equation (^ with the Hamiltonian 

h'o = h{A ~ S)a+a, (13) 
Hint = ^ [(^2 + ^^1 exp {iSt)) a+ -f^ {^2 + exp {-iSt)) a] + hx{a'^af-. 

It is obtained from (^ by the replacements uii(2) ^^2(1)7 ^1(2) ^ ^2(1)7 A ^ A — (5. Nevertheless, these marginal 
cases differ in details, as explained below. 
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In an other scenario of transition to chaos the modulation frequency 5 is varied, with other unchanged parameters. 
In the range of small frequencies (5 <C 7 the modulation of the system is adiabatic. So, in the range of bistability 
the system oscillates between the two possible metastable states. With increasing frequency, at (5 > 7, a strong 
entanglement of these states occurs, and the system comes to chaos. It should be noted that, as we will show below, 
the transition from the regular to chaotic dynamics in the classical system, is marked in the quantum system by an 
increase of the von Neumann entropy, as well as by a strong transformation of the Wigner functions. For the case 
(5^7, the dynamics of the system becomes regular again. 

We give the results of QSD analysis in the regimes of strong anharmonicity, considering the parameters from 
x/7 = 0.7 to x/l = O-l- The former case is strongly quantum mechanical, since the maximum mean number 
of oscillatory excitations is about 10, while the case, x/l = 0.1, corresponds to a quasiclassical regime, when the 
maximum in time oscillatory number is about 130. The truncated basis of Fock number states of harmonic oscillator 
is used for the expansion of state vectors , and an initial vacuum state is chosen. 

A. Time-evolution of the mean oscillatory numbers 

The mean excitation number of double driven anharmonic oscillator versus time interval is of particular interest 
in this paper. Figures 2(a), (b) depict the ensemble averaged mean oscillatory numbers n — (a+a) for both cases of 
regular and chaotic dynamics, which are realized in the classical limit. The parameters for Fig. 2(b) are chosen the 
same as for Fig.l (c). The classical oscillatory number n = \a\^ derived from equation (^ for the same parameters 
as the above chaotic regime, is illustrated in Fig. 2(c). From these figures it is evident that while the classical result 
(Fig. 2(c)) shows usual chaotic behavior, its quantum ensemble counterpart (Fig. 2(b)) has clear regular behavior. 
Thus, the chaotic behavior in the classical model, transforms into the periodic dynamics in the quantum treatment 
which involves ensemble averaging. These results indicate, that quantum dissipative chaotic dynamics is not evident 
on the mean oscillatory number. Now we will show how quantum chaos emerges on both the Wigner function and 
the von Neumann entropy. 

B. Chaos on Wigner functions and quantum interference effect 

We apply the QSD to determine Wigner functions for the quantum states of double driven anharmonic oscillator 
during time evolution. For this we use the well-known expression for the Wigner function in terms of matrix elements 
Pnm = {n-\p\ m) of density operator in the Fock state representation 

W{r,e)=Y,PnmW^n{r.e), (14) 

m.n 

where (r, 9) are the polar coordinates in the complex phase-space plane X = r cos 9^Y = r sin 0, and the coefficients 
Wmn{r, 0) are Fourier transform of matrix elements of the Wigner characterization function pTf 

r -(-l)"A/^e'^"'~"^^(2r)"'-"e-2'-'L™-"(4r2), m > n> 
U(_l)my'^e^('"-")''(2r)"-'"e-2'''L;^r"(4r2), n > mJ 

where are Laguerre polynomials. In our calculation we assume, that the oscillator is initially prepared in a vacuum 
state, and regimes of strong anharmonicity x/l = (0-7 0.1) are realized. 

In Fig. 3 we demonstrate moving of our system from the regular to chaotic dynamics by plotting the Wigner function 
for three values of VI2: ^2/1 = 1 (a), ^2/1 — ^i/l = 10.2 (b), ^2/1 = 20 (c), in a fixed moment of time. The values 
of A/7, x/7j and r2i/7 are chosen to lead to bistability in the model of single driven oscillator {0.2 = 0). 

We can see that for the case of a weak second force Fig. 3(a) the Wigner function has two humps, corresponding 
to the lower and upper level of excitation of anharmonic oscillator in the bistability region. The hump centered close 
to X = Y = Q describes the approximately coherent lower state, while the other hump shows, that the upper state 
is squeezed. The effect of squeezing is displayed as the squeezing of Gaussian. This result represents the known 
property of the Wigner function for the single driven anharmonic oscillator |28j. The graphs in Fig. 3 are given at an 
arbitrary time, exceeding the damping time. As calculations show, for the next time intervals during the period of 
modulation t — 2t:/5, the hump corresponding to the upper level rotates around the central peak. When we increase 
the strength of the second force, the classical system reaches to a chaotic dynamics. The Wigner function for the 
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chaotic dynamics is depicted in Fig. 3(b). Further increasing ^2, the system returns to the regular dynamics. The 
corresponding Wigner function at an arbitrary time exceeding the transient time is presented in Fig. 3(c). It contains 
only one hump, rotating around the centre of the phase-space within the period. As mentioned above, the limit 
^2 3> ill is physically equivalent to the opposite case fii ^ fl2, when the system in each moment of time is close to 
the model of single driven nonlinear oscillator. The difference is that the Rabi frequency fl2 in this case (Fig. 3(c)) is 
taken outside the bistability range, where the system is in the upper level. So, the Wigner function has approximately 
the same form, as the Wigner function for single driven anharmonic oscillator, in the monostable above threshold 
regime of operation, when the system is excited. 

As we see, the Wigner function reflecting chaotic dynamics (Fig. 3(b)), has a complicated structure. Nevertheless, 
it is easy to observe that its contour plots in {X, Y) plane are generally similar to the corresponding classical Poincare 
section. Now we will consider this problem in more detail, comparing the results summarized in Fig.l with the 
numerical calculations of Wigner functions, for the same sets of parameters as for the classical maps. 

We present our results in Fig. 4. It can be seen in Fig. 4. (a) that for the deep quantum regime (x/7 — 0.7, A/7 — —15, 
S/'y — 5), the contour plot of the Wigner function is smooth and concentrated approximately around the attractor 
(Fig. 1(a)). Nevertheless, the different branches of the attractor are hardly resolved in Fig. 4. (a). It can also be seen, 
that in this deep quantum regime, an enlargement of the Wigner function occurs in contrast to the Poincare section. 

Taking a smaller x/7, the contour plot of the Wigner function approaches the classical Poincare section. This can be 
seen in Figs. 4(b), (c). For the last case the correspondence is maximal, and some details of the attractor (Fig. 1(c)) are 
resolved much better in Fig. 4(c). This analysis allows us also to note that the scaling symmetry of strange attractors 
shown on Figs. 1 (a), (c) dissapears for corresponding contour plots of Wigner functions [Figs. 4 (a), (c)] in the 
quantum treatment of dissipative chaos. 

It should be specified that for all contour plots in Fig. 4, the corresponding Wigner functions have regions of negative 
values. These results are related to chaotic regimes. Obviously, this fact reflects a quantum interference in the chaotic 
state of the system considered. It follows that the ranges of negative values of the Wigner function increase with the 
increasing parameter x/7i when the system moves to a deep quantum regime. Nevertheless, the ranges of negative 
values of the Wigner function are also observed for the comparatively small parameter x/7j where an operation regime 
close to the semiclassical is realized. It is illustrated in Fig. 5 for parameter x/7 = 0.1, where the mean excitation 
number equals 130. 

Since our model is dynamical, we can also consider the correspondence between the Wigner function and the 
Poincare section, which are taken at i„ = to + ^y^, for arbitrary tQ. Despite different forms, which the Wigner 
function and the Poincare section have acquired, the correspondence features are the same. 

C. The emergence of chaos on the von Neumann entropy 

As shown in Section IIIA, the quantum chaotic dynamics is not displayed clearly on ensemble averaged oscillatory 
excitation number. The purpose of this subsection is merely to demonstrate how chaos is seen on the von Neumann 
entropy, which is one of the significant characteristics of a quantum ensemble. We also clarify other important 
questions in relation with quantum chaos and entanglement in the system. These analysis complete the above studies 
of quantum dissipative chaos in the Wigner function. 

Two ways of producing chaos in a controlled manner will be considered by monitoring the system through varying 
either the strength of the driving force, or the difference frequency S. We calculate the evolution of the entropy 
by formula (^, using the results for reduced density matrix expressed through an ensemble of trajectories. The 
calculations are performed by diagonalization of the matrix p„,„ in the truncated Fock states basis. 

First, (Fig. 6) we demonstrate the transition of the system from regular dynamics to chaos by plotting the von 
Neumann entropy for three values of strength of the driving field: O2/7 = 1 (a), ^^2/7 = ^^1/7 = 10.2 (b), $^2/7 = 20 
(c). The same parameters are chosen as in Fig. 3 for the Wigner function. In Fig. 6(a) and Fig. 6(c) we plot the entropy 
evolution for the regions with classically regular behavior, while the Fig. 6(b) shows entropy production for the chaotic 
motion. The essential difference between the behavior of the entropy for regular and chaotic dynamics is clearly 
displayed in these figures. The common feature is, that for times exceeding the time scale of transient dynamics, the 
entropy production in the chaotic regime dominates over the entropy production of regular dynamics. These results 
are in good qualitative agreement with the above results on the Wigner function, and reflect the dependence of entropy 
production on both the quantum entanglement, and formation of states in the system. Naturally, entropy production 
is stipulated by the entangling interaction between the anharmonic oscillator and its environment on the one hand, 
and is determined by the structure of mixed states of nonlinear oscillator, on the other hand. As a result, the maximal 
value of the entropy (Fig. 6(b)) is realized for chaotic dynamics with a large number of mixed states as depicted in 
Fig. 3(b), while its minimum takes place for the regular dynamics with one hump Wigner function (Fig. 3(c)). Also, 
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we note, that the oscillations of the entropies evident in Fig. 5, have the frequency 6. The simulations also show the 
definite difference between transient times of regular and chaotic dynamics. There is some ambiguity in the definition 
of the transition time because of the oscillatory nature of curves. Nevertheless, it is clearly evident, that the transient 
time of entropy evolution for the regular case exceeds one for the chaotic dynamics. 

In addition to providing criteria to characterize chaos, we have also studied the behavior of the entropy versus 
controlling parameters Q2 and S. The results of numerical calculations at a definite time moment exceeding the 
transient time, are presented in Fig. 7(a), where minimum values of entropy for different f22 are presented. Let us 
compare these results with the emergence of chaos in the classical limit. It is easy to check, that for the classical 
case and for parameters shown in Fig. 6 chaos appears at the critical point of ^2 = ^cr — 8.195, and disappears at 
il2 — 12.745. It is a well known property of classical chaos, that it appears suddenly. As it is seen in Fig. 7(a), the 
quantum chaos appears smoothly: the entropy increases, as the value of Q2 approaches the critical value ^cr- Another 
way to probe chaos is the varying of detuning S. In Fig. 7(b) the behavior of the entropy versus modulation frequency 
5 are displayed. 

IV. CONCLUSION 

We have presented novel type of a time-dependent systems showing chaotic dynamics and intrinsically quantum 
properties. These systems are modelled by a dissipative anharmonic oscillator driven by two forces of different 
frequencies. We emphasize that this model is different from the ones like single driven nonlinear oscillator, where a 
pulsed pump field could be used and might be proposed as a possible experimental test of quantum chaos in area of 
quantum optics with continuous cw laser. The proposed model seems experimentally feasible with state-of the-art 
equipment and can be realized at least in two experimental schemes. So, the nonlinear behavior of single-mode field 
in a medium with a third-order nonlinearity may provide simple realization of the dynamics of driven anharmonic 
oscillator. In fact, a single-mode field is well described in terms of anharmonic oscillator, and the nonlinear media 
could be an optical fiber or a x'^^ crystal, placed in a cavity. In the later case the anharmonicity of mode dynamics 
comes from the self-phase modulation due to the photon-photon interaction in x'^'^^ -medium, and dissipative effects 
arise from the leakage of photons through the cavity mirrors. Such a system under two driving fields is described by the 
Hamiltonian (|l|) with a, a"*" being the operators of a single cavity mode. Cyclotron oscillations of a single-electron in 
a Penning trap with a magnetic field, is another realization of the quantum anharmonic oscillator. Its anharmonicity 
comes from nonlinear relativistic correction to an electron motion, while dissipative effects arise from the spontaneous 
emission of synchrotron radiation fl^ . The trapped electron driven by a single coherent field, has been experimentally 
realized and studied in Refs. [ p3| . The model we present here corresponds to the one-electron cyclotron oscillator in 
two coherent fields at different frequencies. The corresponding Hamiltonian is given by (|^), where operators a, a+ 
describe the cyclotron quantized motion at the cyclotron frequency. The values of the parameter x/l used in our 
calculations, have been achieved experimentally for both above mentioned physical systems. 

The dynamics of double driven anharmonic oscillator exhibits a rich phase-space structure, including the regimes of 
regular and chaotic motion, with the two Rabi frequencies fli and fl2, and the difference 6 between driving frequencies 
being the control parameters. We suppose that the adequate way of investigating the quantum chaos is not only the 
investigation of the behavior of an individual realization of trajectories, as suggested by several authors, but also 
studying the dynamics of statistical ensemble of quantum trajectories, which is naturally realized in experiments. For 
realization of this program of studies the quantum state diffusion simulation method, based on the master equation of 
Lindblad form is used. We also conclude, that the distinction between regular and chaotic dynamics can be most easily 
understood by studying the dynamics of essentially quantum properties, as the von Neumann entropy and the Wigner 
function are only two examples. In fact, our numerical analysis has shown that the quantum dynamical manifestation 
of chaotic behavior, does not appear on ensemble averaged oscillatory excitation numbers, but is clearly seen on 
the entropy and probability distributions. The connection between quantum and classical treatments of chaos was 
realized by means of comparison between strange attractors on the classical Poincare section, and the contour plots 
of the Wigner functions. We have demonstrated that for small values of the ratio x / 1 the contour plots of Wigner 
functions are relatively close to the strange attractors. Indeed, as we have shown on Fig. 1 (c) and Fig. 4 (c) some 
details of the attractor are resolved on the contour plots. Such likeness of quantum and classical distribution vanishes 
in the deep quantum regime [see Fig. 1 (a), and Fig. 4 (b) for X / 7 = 0.7]. An important point to emphasize is that 
the scaling symmetry of strange attractors in the model studied here violates in the quantum treatment of chaos. The 
drastic difference between the time-dependent behavior of the von Neumann entropy for regular and chaotic regimes is 
clearly displayed in Figs. 6. The von Neumann entropy indicates the connection between chaos and entanglement, and 
has also been used to study characteristic time scales of the emergence of chaos. For short time intervals t < 0.47~^ 
the entropies are linearly increasing function of time for both the regular and chaotic regimes. For times exceeding 
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the time scale of transient dynamics this behavior transforms to the periodic one, however the entropy for the chaotic 
regime dominates over the entropy for the regular dynamics. It has also shown that the transient time for the chaotic 
regime is smaller than in the regular regime. 

In our analysis we have not investigated all possible quantum effects of chaotic dynamics. In particular, we have 
specified that the Wigner function for chaotic regime has regions of negative values even for relatively high values 
of X / 7 (x / 7 — 0-1 on Fig. 5, where the mean oscillatory number n — 130). This fact reflects the quantum 
interference effect in the chaotic dissipative dynamics. However, we have not analyzed the correlation between the 
emergence of chaos and quantum interference, which is an interesting albeit complicated option for the future. 
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FIGURE CAPTIONS 

Fig.l. Poincarc section (approximately 20000 points) for the dimensionless classical complex amplitude of double 
driven anharmonic oscillator, plotted at times of constant phase tnS = 2im. The dimensionless parameters are in the 
region of chaos: 

(a) xh = 0.7, A/7 = -15, I^i/t = f^a/T = 10.2, <5/7 = 5; 

(b) xh = 0.5, A/7 = -25, f^i/7 = ^2/1 = 25, (5/7 = 15; 

(c) xh = 0.1, A/7 = -15, ni/7 = O2/7 = 27, 5/7 = 5. 

Fig. 2. (a), (b) Ensemble averaged (over 1000 trajectories) mean oscillatory numbers (o^a) versus dimensionless 
scaled time for both regular (a) and chaotic (b) dynamics for parameters: (a) xh — O-l, A/7 = —15, ^1/7 = 
1^2/7 = 27, 6/ J = 50; (b) x/7 = 0.1, A/7 = -15, Q^h = ^2/7 = 27, 6h = 5. 

(c) The classical oscillatory mean number jaj for chaotic dynamics. The parameters are the same as for (b). 

Fig. 3. Transition from regular to chaotic dynamics on the Wigner functions for three values of $^2/7 : ^^2/7=1 
(a); ^2/-/ = 10.2 (b); ^2/7 = 20 (c). The parameters are: x/l = 0.7, A/7 = -15, ^1/7 = ^2/7 = 10.2, S/7 = 5. 
The averaging is over 2000 trajectories. 

Fig. 4. Contour plots of Wigner functions corresponding to chaos. The parameters for the cases of (a),(b) and (c) 
are the same as in Figs.l(a),(b) and (c) respectively. The averaging is over 2000 trajectories. 

Fig.5. Illustration of quantum interference effect on the Wigner function for the parameters: x/7 = 0.1, A/7 = —15, 
O1/7 = ^2/7 = 27, 6/7 = 5. The averaging is over 1000 trajectories. 

Fig. 6. Transtition to chaos on the von Neumann entropy: (a) and (c) - regular dynamics, (b) - for chaotic dynamic. 
The averaging is over 2000 trajectories. 

Fig. 7. Behavior of the minimal (during the period) values of von Neumann entropies versus controlling parameters. 

(a) . Dependence of entropy from ^^2/7 for the parameters: x/7 = 0.7, A/7 = —15, ^^1/7 = 10.2. S/7 = 5. 

(b) . Dependence of entropy from S/7 for the parameters: x/7 = 0.7, A/7 = —15, SI1/7 = ^^2/7 = 10.2. The 
averaging is over 2000 trajectories. 
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